Short-range magnetic ordering process for the triangular lattice compound NiGa2S4: a 

positive muon spin rotation and relaxation study 
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We report a study of the triangular lattice Heisenberg magnet NiGa2S4 by the positive muon 
spin rotation and relaxation techniques. We unravel three temperature regimes: (i) below Tc = 
9.2 (2) K a spontaneous static magnetic field at the muon site is observed and the spin dynamics 
is appreciable: the time scale of the modes we probe is ~ 7 ns; (ii) an unconventional stretched 
exponential relaxation function is found for Tc < T < Tcioss where Tcross = 12.6 K, which is a 
signature of a multichannel relaxation for this temperature range; (iii) above Tcross, the relaxation is 
exponential as expected for a conventional compound. The transition at Tc is of the continuous type. 
It occurs at a temperature slightly smaller than the temperature at which the specific heat displays 
a maximum at low temperature. This is reminiscent of the behavior expected for the Berezinskii- 
Kosterlitz-Thouless transition. We argue that these results reflect the presence of topological defects 
above Tc. 

PACS numbers; 75.40.-s, 75.25. +z, 76.75. +i 



On cooling, in the same way as liquids, magnetic ma- 
terials, usually crystallize to form long-range periodic ar- 
rays. However, magnetic materials with antiferromagnet- 
ically coupled spins located on triangular motifs exhibit 
geometrical magnetic frustration which may prevent the 
crystallization to occur [Ij. Such materials are a fertile 
ground for the emergence of novel spin-disordered states 
such as spin liquid or spin glass even without crystalline 
disorder. The simplest example of geometrical frustra- 
tion, stacked two-dimensional triangular lattices with a 
single magnetic ion per unit cell, have been studied in- 
tensively [3]. NiGa2S4 is a rare example of such a two- 
dimensional antiferromagnet which does not exhibit a 
long-range magnetic order at low temperature and which 
is characterized by gapless excitations 0, 0] ■ These phys- 
ical properties lead naturally to the assumption that its 
ground state is a spin liquid. Here we report muon spin 
rotation and relaxation measurements (/iSR) which show 
that a spontaneous static magnetic field appears below 
Tc = 9.2 (2) K where an appreciable spin dynamics is 
measured. In addition, we find the spin dynamics to be 
unconventional in the temperature range Tc < T < Ta-oss 
where Tcmss = 12.6 K. 

Polycrystalline powder of NiGa2S4 has been obtained 
from a solid-vapor reaction using a stoichiometric mix- 
ture of pure elements. The synthesis took place in an 
evacuated silica ampoule. A slow heat treatment over 
several days up to 1000°C has beenperformed in respect 
to the high sulphur vapor pressure [5| . After a final grind- 
ing, the powder was treated for a week at 1000°C. 

NiGa2S4 is a chalcogenide magnetic insulator with the 
Ni^+ magnetic ions (spin S — 1) sitting on an regular 
triangular lattice. The interactions are of the Heisenberg 
type, referring to the isotropic Curie constant fil]. The 



crystal structure consists of two GaS layers and a central 
NiS2 layer stacked along the c-axis. A Rietveld refine- 
ment of a neutron powder diffraction pattern recorded 
at 50 K at the cold neutron powder diffractometer DMC 
of the SINQ facility at the Paul Scherrer Institute (Vil- 
ligen, Switzerland) is consistent with the P3ml space 
group. The lattice parameters are a = 0.3619 nm and c 
— 1.1967 nm, in agreement with previous results fs*]. The 
refinement is consistent with the nominal stoichiometry 
and no impurity phase is detected (detection limit: 1 %). 

Further characterizations of our sample have been done 
by zero-field specific heat and susceptibility measure- 
ments. The magnetic specific heat Cm (divided by the 
temperature) is displayed in Fig. [2l It is in very reason- 
able agreement with the data of Nakatsuji and coworkers 
3]. As these authors reported, we also find that the 
susceptibility recorded under a field of 0.01 T displays a 
weak kink at T^ ~ 8 K. 

Now we report on our zero-field ^SR data; see Refs. 
for an introduction to the /iSR techniques. A spec- 
trum is expressed as aoP^^^{t) where oq is an amplitude 
or asymmetry and P^''^(t) a polarization function. In 
Fig. [1] we display three spectra which probe the three 
temperature regimes that we have unveiled. 

At variance with Ref . j3] , a strongly damped oscillation 
is observed at low temperature where nano-scale mag- 
netic correlations have been detected by neutron diffrac- 
tion [3I; see the spectrum at the bottom of Fig. [TJ Re- 
markably, this oscillation which reflects the presence of 
a spontaneous internal field, vanishes at a temperature 
which is about half the temperature at which the neutron 
magnetic reflections disappear. This is further discussed 
below. Our observation is also technically interesting. It 
shows that, contrary to common wisdom, the detection 
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FIG. 1: (color online). Three zero-field fiSK spectra recorded 
for a powder sample of NiGa2S4 with the GPS spectrometer 
of the S/iS facility at the Paul Scherrer Institute, Villigen, 
Switzerland. The spectra probe three distinct temperature 
regimes. At the bottom we display a spectrum recorded deep 
into the short-range ordered state. The solid line is the result 
of a fit to a function given in the main text. The spectrum 
in the middle has been taken at 11 K where the fit relaxation 
function is a stretched exponential (full line). For reference, 
the dashed line displays the best fit to an exponential relax- 
ation. Such an exponential function is observed at 15 K (top 
spectrum). For clarity, the two high temperature spectra are 
shifted up by 0.05 relative to the third one. Their time scale 
is given by the top horizontal axis, while the scale for the 
low temperature spectrum is at the bottom horizontal axis: 
the range spanned by these two time scales differs by a factor 
~ 30. 



for a material of a zero-field /iSR oscillation is not a fin- 
gerprint of a long-range magnetic order. 

The spectrum at 2.3 K has a steep slope for time 
t < 0.02 IIS. Such a shape is typical for an incommensu- 
rately ordered magnet which shows a characteristic field 
distribution at the muon site 0] . This translates in time- 
space to a Bessel function rather than cosine oscillations 
[6| . Our observation is not surprising since it has been es- 
tablished by neutron diffraction diffraction that the mag- 
netic structure is indeed incommensurate The spec- 
trum is described by the sum of two components for the 
compound and a third component which accounts for the 
muons missing the sample and stopped in its surround- 
ings: ao-Pr^'l^) = aosJo(7/.-Bmaxt)exp(-72A2t2/2) 
arei exp(— A^-^) + abg. Jo is the zeroth-order Bessel func- 
tion of the first kind, i?niax stands for the maximum of the 
spontaneous static local magnetic field distribution at the 
muon site, characterizes the broadening of the probed 
field distribution and 7^ is the muon gyromagnetic ratio 
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The component of am- 



plitude Orel gauges the relaxation of the muons sensing a 
field parallel to their initial polarization. The associated 
spin- lattice relaxation rate is denoted Xz- The asymme- 



try ratio of the first two components is aos/cirei — 1-6 (2) 
to be compared to an expected value of 2. A weak textu- 
ration of the sample could explain this small deviation. 

Figure [2] shows i?max(T')- The full line is a fit to 
the phenomenological function i?max(r) — i?max(0)[l — 
(r/rstat)T with & = 2 and /3 = 0.365. We find 
5max(0) = 222 (20) mT and the temperature for which 
the spontaneous field vanishes Tgtat = 9.3(1) K. The 
magnitude of i?max reflecting the ordered magnetic mo- 
ment, the transition at Tgtat is consistent with a contin- 
uous transition. The field distribution is quite broad- 
ened, A/i?niax — 0.15, in accord with our observation of 
only two oscillations. The detection of these oscillation 
provides a bound for the time scale of the magnetic cor- 
relations [13]: Tc ^ l/(7;iBniax) = 5 US. Last but not 
least, Xz is far from being negligible. Its temperature 
dependence mimics i3max(T') with a value of 12 yits"^ 
at low temperature. Because of the finite Xz value, we 
can derive more than a bound for Tc. From the rela- 
tion Xz = 01, we compute Tc ~ 7 ns, since 
A = 35 (4) mT. This value and the previous bound are 
consistent. Therefore the time scale of the magnetic cor- 
relations is Tc ~ 7 ns at low temperature. The neutron 
scattering results, which set a lower bound of Tmhi = 
0.3 ns for this time, are fully consistent 

The spectrum in the middle of Fig.[T]has been recorded 
at a temperature above Tgtat and in a regime where 
the relaxation is stretched exponential- like: aoP^^{t) = 
Orel exp[— (Azi)'^^] + ^bg- This occurs in the range Tstat < 
T < Tcross = 12.6 K for which three spectra have been 
recorded: at 10, 11 and 11.6 K. We get consistently 
/3^ = 0.56 (3). 

The third spectrum in Fig. [1] corresponds to the case 
where the sample is magnetically homogeneous and the 
magnetic fluctuations are fast enough to give rise to an 
exponential relaxation function i.e. (3^ — 1. As seen in 
Fig. [21 Xz increases for decreasing temperature. This is 
the signature of the slowing down of the magnetic fluc- 
tuations. The relaxation rate is fitted to the formula 
Xz = Xa[T/{T - Tdy„)]2 where Aq = 0.028 (3) fis'^ and 
Tdyn = 9.2 (2) K. This formula is derived from the ran- 
dom phase approximation (RPA) [ll|. It is expected to 
be valid for T > Tgtat provided T is not too close to Tgtat 
which is the case for T > Tcross- As expected, the char- 
acteristic temperature Tdyn introduced in this model is 
equal to Tgtat within error bars: we henceforth replace 
Tdyn and Tstat with a unique parameter Tc — 9.2 (2) K. 

We now discuss the implications of our results for the 
physics of NiGa2S4. 

We note that Cm displays a rounded peak with a max- 
imum at ~ 13 K, i.e. just above Tc- Such a behavior is 
reminiscent of the Berezinskii-Kostcrlitz-Thouless tran- 
sition temperature relevant for a planar XY magnet; see 
Fig. 9.4.3 of Ref. [l2|- This is a key result from our study. 

The relaxation function from a magnetically homoge- 
neous sample is expected to be exponential when mea- 
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FIG. 2: (color online). Upper panel: thermal dependence of 
two parameters extracted from fiSH spectra. On the left-hand 
side is presented the maximum of the static local magnetic 
field at the muon site, -Bmax, for T < Tatat. On the right-hand 
side of the same panel is displayed the relaxation rate, Xz, for 
T > Tcross, i.e. when the relaxation function is exponential. 
The full lines are results of fits described in the main text. 
Lower panel: zero-field magnetic specific heat divided by the 
temperature, Cui/T, deduced from the measured specific heat 
using a dynamic adiabatic technique after subtraction of the 
lattice contribution. The full line is obtained for a linear 
dependence of C'm/T. 



sured in its paramagnetic state; see e.g. elemental Ni [13 1 
or Fe 14 1 or the intermetallics GdNis [l^. A stretched 
exponential relaxation has been observed for a wide va- 
riety of physical quantities in many different systems 
and research areas It arises from a continuous 

sum of exponential decays 
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A square-root relax- 
ation of the Nuclear Magnetic Resonance [l^ and fiSK 
[lit relaxation function is also observed for spatially dis- 
ordered systems. It stems, e.g. for spin-glass materi- 
als, from the distributed nature of the coupling of the 
probe to its environment. Our sample is spatially homo- 
geneous. We ascribe the stretched exponential relaxation 
for Tc < T < Tcross to a multichannel relaxation process. 

Nakatsuji et al have presented the temperature depen- 
dence of the difference AI{T) between the elastic powder 
neutron scattering intensity recorded at T and 50 K, for 
wave vector qc = 0.58 [3|. A/(r) decreases smoothly 
up to ~ 18 K where it becomes negligible. This temper- 
ature is quite different from Tc, where the muon sponta- 
neous field vanishes, and from T^. These differences are 
not surprising, given the different magnetic modes which 



are probed by the three techniques. The neutron inten- 
sity integrates the gc-correlations for times longer than 
Tmin- -Smax IS bullt up from the sum of the dipole mag- 
netic fields at the muon site due to a restricted q range 
of the Ni^"*" magnetic moment Fourier components, and 
characterized by a time scale longer by an order of mag- 
nitude than Tynin- d.c susceptibility probes an even longer 
time scale. That A/(T) vanishes only at ~ 18 K suggests 
a relatively large continuous spectrum of magnetic fluc- 
tuations in NiGa2S4. This could be a common property 
of geometrically frustrated magnets since it has already 
been encountered for Tb2Sn207 [2l[. Because of 

the notable reduction of the Ni^+ {S = 1) magnetic mo- 
ment as measured by neutron diffraction (~ 25%), the 
spectral weight must also extend to time scales smaller 

than Tmin- 

Previous thermodynamic and neutron results [^ and 
the present ^SR measurements on NiGa2S4 demonstrate 
that this two-dimensional Heisenberg triangular lattice 
antiferromagnet has unique low-temperature properties: 
(i) while its Curie-Weiss temperature is rather large, 
©cw = —80 (2) K, it does not display a long-range 
magnetic ordering, but only an incommensurate short- 
range order with a nano-scale correlation length and 
a spontaneous static interstitial magnetic field below 
Tc — 9.2 (2) K; (ii) its ground state is highly degener- 
ate; (iii) there is a slowing down of magnetic fiuctuations 
as the compound is cooled down through Tc rather than 
a spin-freezing as observed for canonical spin glasses. In 
fact, we have found that NiGa2S4 exhibits a conventional 
paramagnetic spin dynamics down to Tcross, where an ef- 
fective multichannel relaxation process sets in down to 
Tc. A finite spin dynamics is detected below Tc. 

Before discussing theoretical proposals for the ground 
state of an exact triangular lattice of isotropic spins in 
light of the experimental results, we note that different 
techniques show that the exchange interaction between 
third neighbors is strong 22, 2^. This result provides a 
clue for the existence of strong frustration, in agreement 
with experimental data [^. 

A hidden order parameter associated to an ordered ne- 
matic phase can be proposed for this system , on the 
ground of the large coherent length inferred from an anal- 
ysis of the specific heat data [3i] . In this case there are no 
long-range two-spin correlations and only the quadrupole 
moments of the Ni^"*" ions order in a long-range manner. 
This phase is classified as a spin liquid. It can be stable 
on a two-dimensional triangular lattice and massless ex- 
citations are present 2^ 2^ 27, 2^. However, at least 
for the model available, strong biquadratic interactions 
are required. 

A second possibility for the ground state relies on the 
presence of residual defects in the system, the effect of 
which is enhanced by frustration. This picture would nat- 
urally explain the stretched exponential muon relaxation 
function for Tc < T < Tcross- Howevcr, the substitution 



4 



of only 1% of Zn for Ni dramatically affects the specific 
heat [J|. Hence, the amount of residual defects is prob- 
ably small. In addition, we are not aware of a model 
calculation which would explain the limited correlation 
length and the persistence of relatively fast fluctuation 
modes below T^. 

A third candidate model attributes the unique proper- 
ties of NiGa2S4 to topological defects inherent to Hcisen- 
berg triangular two-dimensional systems, the so-called 
Z'2-vortices [2^ [3^. Gapless excitation modes and a 
nearly constant susceptibility are predicted ^31,]. Based 
on Monte Carlo simulations, a phase transition was sug- 
gested [291]. However, because of spin- wave interactions, 
the correlation length is finite Jil, jM,] . It is there- 

fore tempting to attribute the transition at Tc to the 
dissociation of the Z2-vortices and Tcross to the crossover 
temperature where the spin dynamics starts to be driven 
by usual Heisenberg spin fluctuations. Remarkably, the 
Z2-vortices manifest themselves in the temperature vicin- 
ity where Cm has a rounded peak. The multichannel 
relaxation for Tc < T < Tcross reflects the magnetic dis- 
order induced by the unbinding of Z2-vortices. 

In conclusion, we have found that, on cooling, the 
frustrated two dimensional triangular lattice compound 
NiGa2S4 first behaves as a conventional magnetic com- 
pound ordering at Tc = 9.2 (2) K. This behavior is 
however observed only down to 12.6 K (i.e. ~ 3.4 K 
above Tc). Below this temperature the fiSH relaxation 
is stretched exponential-like. This result is interpreted 
as the signature of an intrinsic property of the triangular 
system such as the Z2-vortices. The dynamics is never 
frozen, even far below Tc. Finally, the transition at Tc 
from the short range ordered to the paramagnetic phase 
is of the continuous type and occurs at a temperature 
just below that of the specific heat bump, reminding the 
Berezinskii-Kosterlitz-Thouless transition. 

For a further insight into the properties of NiGa2S4, 
it is necessary to examine the wavevector dependence of 
the fluctuating magnetic modes below Tc. On the theo- 
retical front, an interesting result would be to determine 
whether modes with a temperature dependent gap van- 
ishing at T = K are possible for the triangular lattice 
as it seems to be the case for the kagome structure (35| . 
This would provide an explanation for the observed per- 
sistent spin dynamics for T <C Tc (sij . 

Note added. A related report including fiSK data 
recorded on this material is also available [37|. We note 
that the time range available at the /xSR facility used to 
record these data does not allow the authors to evidence 
the spontaneous muon spin precession that we have ob- 
served. 
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